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ABSTRACT 

Context. Debris disc analysis and modelling provide crucial information about the structure and the processes at play in extrasolar 
planetary systems. In binary systems, this issue is more complex because the disc should in addition respond to the companion star's 
perturbations. 

Aims. We explore the dynamical evolution of a collisionally active debris disc for different initial parent body populations, diverse 
binary configurations and optical depths. We focus on the radial extent and size distribution of the disc at a stationary state. 
Methods. We numerically follow the evolution of 10 5 massless small grains, initially produced from a circumprimary disc of parent 
bodies following a size distribution in dN oc s~ 35 ds . Grains are submitted to both stars' gravity as well as radiation pressure. In 
addition, particles are assigned an empirically derived collisional lifetime. 

Results. For all the binary configurations the disc extends far beyond the critical semimajor axis a cri , for orbital stability. This is 
due to the steady production of small grains, placed on eccentric orbits reaching beyond a cri , by radiation pressure. The amount of 
matter beyond a cri , depends on the balance between collisional production and dynamical removal rates: it increases for more massive 
discs as well as for eccentric binaries. Another important effect is that, in the dynamically stable region, the disc is depleted from its 
smallest grains. Both results could lead to observable signatures. 

Conclusions. We have shown that a companion star can never fully truncate a collisionally active disc. For eccentric companions, 
grains in the unstable regions can significantly contribute to the thermal emission in the mid-IR. Discs with sharp outer edges, 
especially bright ones such as HR4796A, are probably shaped by other mechanisms. 
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1. Introduction 

1.1. Planets and discs in binaries 

Debris discs are dusty circumstellar systems observed around 
main sequence stars. They are thought to represent a late stage 
of disc evolution, when the primordial gas has been dispersed 
and most of the solids have been accre ted onto embryos or fully 
formed planets (Lagra nge et al.l l2000). What is seen are small, 
typically < 1 cm dust particles, produced by collisions among 
larger unseen parent bodies, probably leftover pla netesimals no t 
used in the planet building process (see review bv lWvatt 1 ^008). 
Most debris discs are initially detected through an excess IR 
flux in the stellar spectrum due to the thermal emission of dust. 
In some favourable cases resolved images are also obtained, in 
scattered light and, for some systems, also in thermal light at 
longer wavelengths. Most resolved discs do not appear homo- 
geneous and smooth but show on the contrary pronounced spa- 
tial features, l ike clumps, warps, radial asymetries or spiral arms 
( Wvatt 2008). Numerical and analytical models have been used 
to analyse these features and have provided precious information 
about the structure and the dynamics of these systems. In many 
cases, the presence of azimuthal, radial or vertical structures has 
been successfully explained as been due to the gravitaional pull 
of unseen perturbers, such as planets or planetary embryos. Q 
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1 The best known examples of model-infered pertu rbers are the 
Jovian planet producing the warp in the /JPictoris disc ( Mouillet et al. 



In this respect, binaries are of special interest, since the sec- 
ondary star's perturbations should significantly affect the ex- 
tent, structure and evolution of a circumprimary debris disc and 
could induce pronounced features: confined rings, spiral arms, 
etc. More generally, studying dusty discs in binaries, especially 
young and bright discs that should be the observable counter- 
parts of the last stages of planet formation ( Ken von & Brom ley 
2004), can provide crucial information about planet formation 
processes in these systems. Even fainter discs around older stars 
are of great interest since they should trace the location of as- 
teroid or Kuiper belt-like belts, whose structure should provide 
clues to the presence of unseen planets. Comparing debris disc 
structures in binaries to that of single star systems could thus 
provide indirect but crucial information about how binarity af- 
fects planet formation processes. 

Such information is very timely in the current context of a 
renewed interest for the problem of planet formation in binaries. 
Recent years have indeed seen a spectacular surge in the num- 
ber of studies dedicated to this issue, fueled by the discovery 
of ~ 80 exoplanets in binaries, some of th em in tight systems 
having a separation of less than 30 AU dDesidera & Barbieril 
120071: Mugrau er & Neuhausen 1 l2009h . Some issues are today 
relatively well understood, for instance the limiting distance 
a cr i,, as a function of the binary 's configuration, for or- 
bital stability around each star (e.g iHolman & Wiegertl [l999: 



Il997h and the one produci ng the inner edge of Fomalhaut's magnifi- 
cent ringjjKalas et al. 2005) . Both planets have since then been directly 
imag ed ( Kalas et al. 2008 ; Lagra nge et alll2009h . 
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iMudrvk & Wul 120061) . This limit does also roughly corre- 
spond to the region where the last stages of planet forma- 
tion, leading fro m lunar-sized embryo s to fully formed plan- 
ets, can proceed (Ouint ana et al.l [2007h . However, some issues 
are still actively debated, in particular regarding the earlier 
stage lead ing from kilom etre-sized planetesimals to the em- 
bryos (e.g.lThebault et al.ll200d 120091: iPaardekooper et al.l2008t 
IXie & Zhoull2009HPaardekooper & Leinhardtll2010l) . 

As signposts of planet formation and/or planet location, de- 
bris discs are thus of crucial importance to understand what is 
going on in potentially planet forming binaries. However, to un- 
derstand the connection between an observed dust population 
and the underlying architecture of the whole disc (collisional 
parent bodies, planets, etc.), it is crucial to first understand what 
imprint the binarity itself will leave on the disc. 

1.2. Previous studies 

Photometric surveys, looking at mid-IR excesses of debris discs 
around large stellar samples, have been attempting to identify 
trends specific to binar y systems. Th e existe nce of such trends is 
however still debated. iTrilling et al.l (I2007I) looked with Spitzer 
at 69 A- and F-type double stars and found that the disc fre- 
quency is lower for binaries with intermediate separation 3- 
30AU than for both ti ghter and wide r binar ies. On the contrary, 
iPlavchan etail d2009l) and iDuchend (|2010|) claim that no such 
trend is visible. Clearly, additional observations are needed to 
settle this issue. In any case, the theoretical part of Trilling et al. 
d2007l) 's study is of great interest. It argues that, since thermal 
emission at a given wavelength A j probes a given dust tempera- 
ture and thus a distance a,/, from a star of a given spectral type, 
then binaries having their a crlt < a t f, < a c f t should have a lower 
excess rate at At, because regions that are probed correspond to 
unstable orbits (where a c % is the inner limit for stable circumbi- 
nary orbits). In any case, the fact that binaries do have an effect 
on debris discs would seem a logical result in the light of re- 
cent studies that have shown that binaries affect the evolution of 
primo rdial protoplanetary discs in younger systems dCieza et al.l 
I2009T) . 

Regarding spatially resolved systems, the most famous ex- 
ample of a debris disc in a binary is probably HR4796A, with 
a very bright ring at ~ 70 AU from the central star and a phys- 
ically bound co mpanion at a projected distance of ~ 510AU 
dJura et alj|1993l) . The possible impact of the companion on the 
disc has often been mentioned as a potentiall y important ef- 
fect without having been inv estigated in details ( Augereau et alJ 
1 19991 ISchneideretaill2009l) . So far, the only system for which 
the effect of a companion star has been numerically investigated 
is HD 141569A. This AOV star is surrounded by a massive dusty 
and gaseous disk with a complex structure, including a spiral 
arm (IClampin et al.ll2003l) . It has two possibly bound M dwarf 
companions at projected distances of 750 and 890 AU, whose 
infl uence on the dust disc ha s bee n studied in a couple of stud- 
ies. lAugereau & Papaloizoul d2004l) have estimated the steady 
state reached by the disc after repeated passages of a compan- 
ion of varying eccentricity eb and separation a/,, but neglect- 
ing the crucial effect of radiation pressure (see Sec ll.3l ). The 
same collisionl ess and radiation-p ressure-less problem was in- 
vestigat ed by lOuillen et alJ Q2005) using a hydrodynamical ap- 
proach. lArdila et aH d2005h have considered the response of a 
dusty disc taking into account radiation pre ssure, but only for a 
single fly-by of a companion. More recently. iReche et al] d2009) 
also looked at the disc's response to a single fly-by. 



1.3. Present approach 

As can be seen, the debris-disc-in-a-binary problem has so far 
not been thoroughly investigated, i.e., only for one speficic sys- 
tem and each time only considering one specific aspect of the 
problem. We intend to take these studies a step further and con- 
sider this problem from a general perspective, exploring a large 
range of disc profiles as well as binary parameters (separation a/, 
and eccentricity eb). We will especially focus on one crucial ef- 
fect, that of radiation pressure on the smaller grains. This mech- 
anism might indeed greatly affect the response of the system to 
a companion star's perturbation. We can basically identify two 
potential effects: 

- 1) Small grains are steadily produced by collisions from 
larger parent bodies, and are automatically placed on high- 
eccentricity orbits by radiation pressure. So that even if the 
parent bodies stay within the orbital stability region, high-/? 
grains might populate dynamically "forbidden" regions, de- 
pending on how fast they are produced as compared to how 
fast they are dynamically ejected. 

- 2) Because most small grains are on high-e orbits with their 
apoastron potentially in unstable regions, they might be re- 
moved by the companion star even if they have been pro- 
duced from parent bodies on stable orbits. In other words, 
companion star's perturbations could remove small grains 
from regions within the stability limit. 

The purpose of this study is to quantitatively investigate these 
different effects. We use a specifically developed numerical 
model to follow the evolution of a collisional debris disc, orbit- 
ing a primary star and perturbed by a more distant companion. 
Grains are submitted to the radiation pressure of the primary star. 

2. Model 

2.1. Methodology 

The system we intend to investigate is that of a debris disc or- 
biting a primary star and perturbed by a stellar companion. The 
dynamical evolution of the debris disc will be traced by a popula- 
tion of N test particles, using a A' 1 ' order Runge-Kutta integrator 
with adaptative step. We adopt the restricted 3 body problem, 
where particles are massless and only the potential of each star 
is considered. One crucial additional force is taken into account, 
i.e., radiation pressure from the central star. With such a deter- 
ministic approach it is impossible to realistically follow the col- 
lisional evolution of the system, since the numerous fragments 
produced by each individual impact would lead to an exponen- 
tial increase of N rapidly exceeding any manageable numberQ. 
However, collisions cannot be fully ignored for the present prob- 
lem. A first reason is that they are imposing the size-distribution 
of the debris population. This parameter is here crucial since 
radiation pressure is a size-dependent force, whose global ef- 
fect is related to the fraction of small, radiation-pressure affected 
grains. In addition, collisions set the rate at which grains are pro- 
duced, a parameter that is in direct competition with the rate at 
which they can be removed by dynamical effects. 

In practice, we first run a simulation with Npb = 10 5 "parent 
bodies" not affected by radiation pressure (J3 = 0) and wait until 
it reaches a steady state. By steady state, we mean that the spatial 



2 a realistic treatment of collision outcomes would require a statis- 
tical approach, wi t h no or a simplified treatmen t of the dynamics (e.g. 
iKrivov et alJl2006l : lThebault & Augereaull2007h 
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structure of the parent body disc is the same between two con- 
secutive passages of the companion star at the same phase angle 
(the disc structure is the same at a time t and at time t + tbm). We 
then consider a set of 10 parent body disc configurations at 10 
regularly spaced time intervals wit hin one bina r y orbi tal period. 
Following a procedure similar to lArdila et"al] d2005l) . we gen- 
erate, from each parent body configuration, a dust population 
follow ing a collisional equilibrium size distribution ( Doh nanvil 
119691) in dN oc s - i 5 ds (i.e., dN oc p L5 d0), from p = 0.012 
down to very small unbound grains with p = 2.5, and we com- 
pute its evolution in time Q. Dust particles are progressively re- 
moved either by radiation pressure for grains with p > 0.5, or by 
collisional destruction or by dynamical perturbations from the 
companion star. For the dynamical ejection, we remove particles 
when they pass at less than 0. 1 Hill radius of the companion, or 
when they are on an unbound orbit and at a radial distance ex- 
ceeding 50a cr ; f . As mentioned earlier, collisions cannot be prop- 
erly taken into account with our deterministic code, so we con- 
sider a simplified approach in which, during each timestep dt 
spent within the collisionally active parent body region, each 
particle is assigned a collisional destruction probability 



Table 1. Set up 
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where tocoll 1S a characteristic collisional destruction timescale. 
We assume that most collisions are destructive, so that tocoli 1S 
comparable to the coll i sion t imescale t c „u for which we follow 
Thebault & Augereau 1 12005) and assume 
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where so is a reference size, 6V and 6Vo are the departures from 
the local Keplerian velocity Vxep of the considered particle and 
for a parent body with p = respectively, and f co // = t or b/(4-n t) 
is the classical size-independent simplified expression, propor- 
tional to the local orbital period t or b divided by the local total 
cross sectional area, or vertical optical depth r. Q. For so we 
consider the size of the smallest particles that behave like par- 
ent bodies, assumed to be the size for which the eccentricity 
imposed by radiation pressure ep = /?/( 1 . - p) is smaller than 
the eccentricity assumed for the parent bodies eo- For the size 
dependence index a, we assume the value -0.3, as obtained by 
IWvatt & Dent! (120021) . Note that the collision rate is assumed to 
be independent of the companion star's orbit and of the dynam- 
ical perturbations it triggers in the birth ring. This simplification 
is acceptable as long as the particle-in-a-box approximation for 
estimating encounter rates holds. However, this approximatio n 
breaks down for particles on high-e orbits (IWvatt et al.l 12010). 
which could be the case for the highest eh considered here, i.e. 
0.75, for which orbits in the parent body disc might reach values 
as high as ~ 0.3. Such values correspond to the li mit where a 
more refined model like that of IWvatt et alJ d2010h is required. 
Although implementing such a model is beyond the scope of the 
present work, we would thus like to point to possible limitations 



3 Note that the size distribution in real systems with a minimum size 
cut-off should depart from this power-law and exhibit a characteris- 
tic wa vy profile (e.g Campo Bagatin et al. 1994; Thebault & Augereau I 
120071) . However, we shall ignore these refinements here and simply as- 
sume the Dohnanyi profile 

4 The apparent difference with IThebault & Augereau! d2005l) 's ex- 
pression comes from the fact that the t co u n considered here implicitly 
regroups several parameters of Thebault & Augereau ( 2005)'s Equ.6 



Radial extent^ 0.6 < a < 1 (extended disc case) 

0.55 < a < 0.6 (inner ring case) 

Eccentricity of the parent bodies e 0.05 

Surface density profile (parent bodies) X = est 

Disc optical depth for t cM estimation r = 2 x 10~ 3 (dense disc) 

r = 2 x 10~ 4 (faint disc) 

Number of test particles N 10 5 

Size ranged /?( s max ) = 0.012 < p(s) < p(s, m „) = 2 

Size distribution at t = dN(s) oc s~ 3S ds 

" normalized to a,„, 



as parameterized by /3 oc \ js 



for our t co u estimate in the high e B cases (for more discussion on 
the collision rate approximation, see also Sec l3.3l >. 

We then follo w the sa me procedu re as in 

IThebault & Augereau I d2005l) and iKrivov et alJ §009) and 
record, at regularly spaced time intervals, the positions of all 
remaining particles. The disc's spatial structure is progressively 
obtained by adding up these instantaneous snapshots until all 
particles have been removed either dynamically or collisionally. 
This procedure is implicitly equivalent to assuming an arbitrarily 
co nstant dust production rate f rom parent bodies (see discussion 
in IThebault & Augereau 1 120051) . By this method we obtain 10 
disc profiles corresponding to our sample of 10 different initial 
parent body configurations. The final disc profile is obtained by 
averaging these different profiles. 



2.2. Numerical set up 

For the sake of clarity and so that the results can be easily under- 
stood in terms of departures from the purely gravitational case, 
all distances will be scaled to a e ,.„, the estimated outer limit for 
orbital stability corresponding to a given binary configuration 
ab, eb,p in the purely gravitational case (i.e., when negecting ra- 
diation pressure, see Sec|l.l|>. We compute q cnr using the empir- 
ical formulae bv lHolman & Wiegertl (1 19991) 

= 0.464-0.380/z-0.631e fc +0.586jue 6 -i-0.150e?-0.198jue?(3) 

ab 

where yu = m.2l{m\ + mj), with m\ the mass of the primary and 
;«2 that of the companion star. Note that there exists other esti- 
mates of this parameter, but the iHolman & Wiegertl (1 19991) es- 
timate has proven to be relatively accurate for most cases and 
is thus satisfying as a first order approximation of the width of 
the dynamically stable circumprimary region. In any case, let us 
stress that our results do not depend on the assumed value for 
a cr u, since we automatically obtain our stable parent body disc 
from specific test runs Q. The a crlt parameter is just used as a 
convenient, and dynamically relevant scaling factor. 

Given the number of free parameters, a thorough phase space 
exploration is impossible. We cannot in particular explore all 
possible disc and binary configurations. We will thus restrict 
ourselves to two disc geometries: 1) A nominal "wide disc" 
case where the parent body extends up to a e ,„, 2) an "inner 
ring" case, whose outer edge lies far within the orbital stability 
limit. For the disc's optical depth, we also consider two cases: 
1) a nominal case with a dense and collisionally active disc 



the accuracy of the lHolman & Wiegertl (l 999) estimate is confirmed 
by the fact that, for all our parent body runs, the maximal extent of the 
parent body disc is in good agreement with a crit 
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Fig. 1. Run with no companion star. Normalized and azimutally 
averaged radial profile of the vertical optical depth for the whole 
grain population (solid line), and for three sub-populations cor- 
responding to different grain size (i.e., /3) range. The initial par- 
ent body disc (J3 = particles) extends to a en , = 1 . 

with t — 2 x 1CT 3 corresponding to bright debris discs such as 
HD32297 or (3 Pictoris (the brightest systems, such as HR4796 
have even higher t, see Sec J4. 1 .2b . and 2) a low r = 2 x 1CT 4 
case modelling fainter systems, with lower collision rates, such 
as AUMic. For each disc case, a full set of binary configuration 
cases will be run. We explore a wide range of e/, values, ranging 
from to 0.75. For each et,, the value of at, is then automati- 
cally given by the condition that a crlt is kept constant. The only 
parameter that will not be thoroughly explored is the mass ra- 
tio, for which we assume a fixed value fi = 0.25 (except for the 
specific HR4796A run, see Sec l4.1.2l >. which corresponds to es- 
timates of the mean mass ratio in binaries derived from exte nsive 
surveys dKroupa et al.lll990l : lDuquennov and MavorHl99Tl) . 

For the dust population, we consider N = 10 s test particles, 
launched from the parent body ring following a dN oc s~ 3 5 ds 
size distribution. The particles' absolute sizes are of no impor- 
tance here, the crucial parameter being their /3 value. As a con- 
sequence, sizes will be scaled by this f3 parameter, and we take 
Pmin = 0.012 (maximum size) and /3 max - 2.5. Due to the steep 
size distribution, we run in practice 3 separate runs with N - 1Q 5 
particles, for the 0.0012 < j3 < 0.083, 0.083 < p < 0.5 and 
0.083 < P < 2.5 domains respectively, which are then weighted 
and recombined. 

All initial conditions are summarized in TabQ] 

3. Results 

3.1. Test run for single star 

We first run a test simulation without a binary companion, i.e., 
just consisting of an initial disc of parent bodies extending to 
o-crit = 1 (flcrit is here only a reference for comparison with the 
cases with companion where it has a physical meaning). Its only 
evolution is thus due to the analytical collisional removal pre- 
scription. FigQ] shows the steady state reached by the system. 
Not surprisingly, the region beyond a crl , is populated by small, 
high-/? grains launched on high eccentricity orbits because of 
radiation pressure. The optical depth profile beyond the parent 
body region has a radial dependence that follows a power law 



Fig. 2. Azimutally averaged radial profile of the vertical optical 
depth for 6 different cases. All profiles have been normalized to 
their value at a — a cr i,.) 
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in oc r~ . This is a well known result for collisional systems at 
steady state, for which it has been shown that the vertical optical 
depth naturally assumes a r~ L5 profile (r -3 - 5 for the mid-plane 
lumi nosity) beyond the coHis i onaily active region of parent bod- 
ies (IStrubbe & Chiand I2006I: iThebault & Wull2008h . Note that 
the contribution of unbound grains (J3 > 0.5) is negligible in the 
whole considered radial region (see Sec l3.3l a more detailed dis- 
cussion). 



3.2. Perturbed disc structure in the presence of a companion 
star 

In the following simulations we include the companion star and 
explore several binary configurations, by varying et, and at, while 
keeping a crll — 1 as a reference anchor point. The optical depth 
t for the disc is for most cases that of our nominal dense disc 
(high collision rate) case, but the faint disc t — 2 x 10~ 4 case 
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Fig. 4. Normalized synthetic images in scattered light, at steady state, for a wide disc of parent bodies initially extending to «„•„, 
and for different orbital configuration of the binary, where the value a c „, is kept constant and <?/, varies. The upper left panel is a 
test case with no companion and the next 5 runs are for =0, 0.1, 0.3, 0.5 and 0.75 respectively. The location of the a crU radius is 
shown by the dashed black circle. The binary configuration is schematized in the small white square at upper right of each panel, 
where the region within a cn/ is shown as a filled black circle. All distances have been normalized so that a en , = 1. 
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is also explored. The main outcomes are summarized in Fig|2] 
which displays the averaged radial profile of the vertical optical 
depth, at steady state, for all cases. 

The main result is that, for all cases, the dusty disc always 
extends well beyond the purely dynamical orbital stability ra- 
dius a cr i t . This is due to the steady production of small, high-/? 
grains that are produced within a crn but are launched on eccen- 
tric orbits reaching far beyond this limit. The dust population in 
the dynamically unstable region is larger and extends farther out 
for higher eccentricities of the binary. In effect, for a fixed a crlt , a 
companion on a circular orbit is the most effective in truncating 
the circumprimary disc, while eccentric companions have a more 
limited influence. This is because perturbers on eccentric orbits 
spend a large fraction close to their apoastron far away from the 
primary, thus allowing small grains that are produced in the cir- 
cumprimay disc to move towards the outer parts of their (radia- 
tion pressure induced) eccentric orbits without encountering the 
secondary star. True, most of these grains will eventually be re- 
moved by the companion's perturbations, but this removal will 
take more time than for a circular orbit companion that always 
moves in the vicinity of the primary. This increased removal time 
allows for more grains to be present at a given time, since the 
collisional production time is, to a first order, independent of the 
binary configuration. In addition, the contribution of unbound 
grains (J3 > 0.5), which is negligible in the no-companion case, 
becomes significant in the r > a cr u regions, as can clearly be seen 
in Fig[3] This is because these grains, contrary to all the others, 
do not get further depleted by the presence of the companion star 
(see also Sec l3.3l ). 

Of course, even for the highest eb cases, there is always 
less dust beyond «„-,, than for the companion free case, as is 
clearly illustrated by Fig|2] showing the averaged radial pro- 
file of the vertical optical depth. Taking as the reference the 
ratio T(2fl m ,)/T( flmI ), we see that, for the et, = case, it can be 
a factor 10 lower than for the companion free profile where 
T(2a cH ,)/T(a cri ,) - 0.1. For the more eccentric binary cases, how- 
ever, the difference with the no-companion case is much more 
limited. It is only a factor 3 below the companion-free case 
for the et, = 0.5 run and less than a factor 1.5 for e/, = 0.75. 
Note that, for all cases, the optical depth profile within the par- 
ent body disc slightly increases from 0.6a cn , to a crn . This is due 
to the contribution of particles with high /? that are more nu- 
merous at the outer edge of the ring, where a large fraction of 
high-/? grains produced further in can be observed, than at the in- 
ner edge, where only locally produced high-/? grains contribute. 
When only considering the contribution of large particles with 
low /?, we observe a flat optical depth profile similar to that for 
the input parent bodies (see FigsQ]and[3]l. A logical consequence 
is that profiles in the birth ring for efficiently truncated discs 
(those with low e/, companions) are flatter than for less efficently 
truncated ones, because in the former cases the birth ring is more 
strongly depleted of small grains (and thus more dominated by 
low /? particles) than in the latter ones (see Sec l3.3l >. 

For illustrative purposes, we show in Fig|4] synthetic head- 
on seen images, in scattered light assuming gray optical scatter- 
ing for the grains, for a set of runs with e/, ranging from to 
0.75, as well as for the no-companion case for comparison. Note 
that these "images" cannot reproduce fine spatial structures, such 
as spiral arms or sharp azimuthal asymmetries, because of the 
phase averaging procedure used in producing the steady-state 
distribution. However, they clearly show that resolved observa- 
tions with a sensitivity range of 10 2 could reveal flux coming 
from dynamically forbidden regions. 
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Fig. 5. Differential fi(s) (i.e., size) distribution within the parent 
body disc for the no-companion and the eb - 0.3 companion 
cases, as well as for the same e/, = 0.3 case with a lower collision 
rate (t = 2 x 10~ 4 ) 



3.3. Size distribution in the parent body region 

Fig|5] shows the steady state size distribution within the par- 
ent body region for both the non-companion and the reference 
eb = 0.3 cases. For the companion-free case, the differential size 
distribution has some small wiggles in the 0.2 < /? < 0.5 domain 
but roughly follows a power law in dN oc /3 l 6 df3, correspond- 
ing to dN oc s~ 36 ds, which is relatively close to the input size 
distribution. More importantly, we do not observe a depletion of 
small, high-/? grains within the birth ring, despite of the fact that 
these small grains are on high eccentricity orbits and spend most 
of their time outs i de of the ring. This confirms th e results of 
IStrubbe & Chiang! d2006h and lThebault & Wul d2008h that debris 
discs tend toward a state where the collisional equilibrium distri- 
bution (be it in s -3 - 5 or not) holds within the collisionally active 
region, meaning that for the system integrated distribution, small 
high-/? grains are in excess. Another expected result is the sharp 
transition at the blow-out size (J3 = 0.5), with a strong depletion 
of unbound grains. This is because these grains leave the system 
on a dynamical timescale while bound particles have a lifetime 
imposed by collisional effects. For the considered value of the 
system's optical depth, a few 10~ 3 , the former is much shorter 
than the latter. 

The effect of the companion on the size distribution is to de- 
plete it from its smallest particles. For the eb = 0.3 case, this de- 
pletion becomes significant for/? > 0.15 grains, while /? > 0.25 
particles are almost completely absent. Note that this depletion 
occurs in the parent body disc, within the orbital stability limit 
a C rit- It is due to the fact that small particles produced at r < a crlt 
are placed by radiation pressure on eccentric orbits that enter the 
potentially unstable region beyond that limit. Depending on their 
position relative to the companion while in the r < a cra region, 
some of these small grains will be ejected and will not come 
back at their periastron in the parent body ring. As could be log- 
ically expected, this small grain depletion gets more pronounced 
when decreasing the disc's optical depth (t = 2 x 10~ 4 case), 
because the collision production rate of small grains is lower, 
while the removal timescale, only due to the dynamics, remains 
the same. Another noteworthy result is that, for the same given 
optical depth, the number of unbound grains is approximately 
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the same with or without the secondary star. This is a logical re- 
sult in itself, because the removal rate of these small grains, im- 
posed by radiation pressure blow out, is only marginally affected 
by the presence of the companion. An interesting consequence 
is that there is no longer a sharp transition at f3 = 0.5. The ex- 
cess of grains just below the j3 = 0.5 threshold has been erased 
by the fact that bound particles in the 0.15 < P < 0.5 range are 
strongly depleted by the secondary star. This is why the contri- 
bution of unbound grains to the total optical depth becomes non 
negligible, as has been seen in Figj2] 

Note that we do not observe here the "wavy" size dis- 
tribution, starting at the blow-out size (J3 ~ 0.5) imposed 
by radiation pressure, expect ed in realistic collisional systems 
( Theb ault & Augereau 1 120071. This is because our expression 
for t co u is too simplified to account for such complex effects. In 
order to do so, we c ould in principle use the empir ical expression 
derived in Equ.7 of Theb ault & Augereau I d2007i) . Nevertheless, 
this expression is only valid for isolated debris discs where the 
collisional cascade is not affected by dynamical perturbers 0. 
This is clearly not the case here, because the companion star 
strongly modifies the abundances of grains close to the blow- 
out size. It tends in particular to erase the discontinuity at the 
f3 = 0.5 limit that is the very source of the wavy distribu- 
tion. We would thus expect a much less pronounced oscillation 
in the size distribution, although a secondary wave could pos- 
sibly start at the limiting size for the depletion by the binary 
(J3 ~ 0.15). Studying these effects would require to use a code 
that can self-consistently follow the coupled dynamical and col- 
lisional evolution of the system. Such a task exceeds by far the 
scope of this paper, as there is to our knowledge no such all- 
encompassing code available yet. FigO should thus be regarded 
as a first approximation. However, to check the robustness of 
our results we performed a te st run using the t co u expression of 
Theb ault & Augereau I (120071) for the — 0.3 companion case. 
We verify that, even for this highly unrealistic case artificially 
forcing an effect at the /3 = 0.5 limit, the resulting size distribu- 
tion is remarkably similar to that of Fig|5] The main part of the 
wavyness being completely erased by the dynamical depletion 
in the f3 ~ 0.15 domain. We are thus relatively confident that our 
main result, i.e. the binary-induced cut-off at /3 ~ 0. 15, should be 
relatively robust given its amplitude, which exceeds by far that 
of any wa vyness observed in detailed collisional studies (see for 
instance iKrivov et al.ll2006t iThebault & Augereau Il2007h . 

3.4. Inner ring case 

We consider here a parent body ring extending only to a OM = 
0.6a enr . This corresponds in practice to astrophysical cases with 
greater binary separation, where the companion star didn't trun- 
cate the initial protoplanetary ring and a ou , is the natural disc 
outer limit. As it can be seen in Fig|6] even in this case, the debris 
disc extends well beyond a eri , for all explored cases. However, 
the amount of material in the unstable zone is smaller than in 
the case where the disc extends all the way to a cr i t . Taking again 
as a reference the ratio T(2a crlt )/T(a crlt ), we see that it is one order 
of magnitude lower than in the unperturbed case (no companion 
star) for binaries with <?/, < 0.5. For systems with higher eccen- 
tricites, however, the difference with the binary-free case drops 
to less than 50%, a value comparable to that of the previous case 
with an extended disc. As for the size distribution within the par- 
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6 When using the expression of Thebaul t & Augereau 1 1[2007) instead 
of our Equ[2]in the non perturbed single-star case we do indeed obtain 
a pronounced wavy pattern for the size distribution 
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ent body ring (Fig|7j), results are remarkably similar to those with 
ciout - acrit- It is depleted from most high-/? grains, the limiting [3 
value being only slightly higher, ~ 0.2 instead of 0.15. 



4. Discussion 

4. 1 . Dust in dynamically unstable regions 

One important result of our investigation is that a companion 
star is never able to fully truncate a collisionally active debris 
disc down to the critical semimajor axis for dynamical stability. 
More specifically, there is always dust in the dynamically unsta- 
ble regions of the system, mostly small grains placed on bound 
eccentric orbits by radiation pressure. Exactly how much mate- 
rial depends on the binary configuration as well as on the disc's 
optical depth. A faint, less collisionally active disc is less able 
to fill the r > a crll region than a denser, more massive system 
(see Figf2]i. Note however that most resolved debris discs have 
optical depths close to, or even higher than our dense disc case 
with t = 2 x 10~ 3 . 
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Fig. 8. Synthetic Spectral Energy Distribution (SED) derived for the extended disc case, for two different spectral types for the 
primary: F5V (upper panel) and A3V (lower panel). For each case, two different <?/, values, 0.3 and 0.75, are considered. The 
companion's semi-major axis is given by the condition that the orbital stability limit is at a cr u = 4 AU, which leads to at, — 16 AU 
for the ei, = 0.3 run and a/, = 50 AU for the et, = 0.75 one. The respective contributions from r < a crlt and r > a cr i t ("unstable") 
grains are plotted. For the A3 case, we also plot the SED expected for a reference unperturbed annulus located at r — a cra . 



4.1 .1 . Signature in thermal emission 

At first g lance, this result cou ld help explaining a puzzling re- 
sult from Trilling et al. (2007), who detected, for at least 3 bi- 
naries, dust in dynamically "forbidden" regions. However, the 
issue might be more complex than that. The first problem is 
that, as acknowledged by the authors themselves, for all 3 cases, 
flux excess is only detected at 70pm, so that only a maximum 
dust temperature Tj can be derived (taking the non detection at 
24pm as a lower limit), and thus a minimum distance to the star. 
There is thus a possibility that the observed dust is further out, 
possibly on a stable circumbinary orbit (see Fig. 5 of that pa- 
per). Furthermore, the non-detection at 24yum should preclude 
the possibility that there is a massive disc closer to the primary, 
which should in principle be the case because the detected "un- 
stable" dust has to be produced from collisions from more mas- 
sive parent bodies further in. On a related note, even if multi- 
wavelengths observations were available, they might not yield a 
characteristic Tj corresponding to dust beyond a enr , because the 
total thermal flux might remain dominated by grains in the stable 
parent body region where the optical depth peaks. 



We investigate this issue by deriving synthetic Spectral 
Energy Distributions (SED) from the steady-state optical depth 
distributions obtained in our dynamical runs. To this effect we 
need to specify parameters that did not need to be quantified in 
our numerical study. The first one is the absolute spatial scale, 
because we need real distances to derive grain temperatures and 
thermal emission. We consider here a fixed a cr u =4 AU value, 
and explore two cases for the binary eccentricity, <?/, = 0.3 and 
et, — 0.75. This gives a binary separation of 16 AU in the first 
case and 50 AU in the second one. These separatio ns are typical 
for the intermediate separations binaries for which Trillin g et alj 
(2007) detected grains in the "forbidden" regions. We assume 
hard sphere silicate grains and calculate their wavelength- and 
size-dependent optical prope rties with the Mie theory (optical 
constants from lDraindl2003l ). We consider two possible stellar 
types, F5 V and A3 V, for which the stellar photospheres are mod- 
elled with synth etic NextGen stellar atm osphere spectra for rep- 
resentative stars (lHauschildt et al.ll999h . At any wavelength, the 
total disc flux, assuming it is unresolved, is obtained by sum- 
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ming up over the distances the thermal emission of all grain sizes 
weighted by the size distribution. 

For an F5 star (Fig[8]A and B), we see that the contribution 
of unstable grains, although non-negligible (up to 15% of the 
flux in the et = 0.3 case and 30% in the et = 0.75 one), is not 
enough to drastically modify the total SED. This is mainly be- 
cause grains in the unstable regions, which have a size < 3 s« ow , 
are too small (stiow ~ l/" m f° r an F5V star) to significantly con- 
tribute at wavelengths > 20/mi where, for an F5 central star, 
their thermal emission should peak if they were perfect emitters. 
For the A3 case (Figj8]C and D), on the contrary, the effect of 
unstable grains is much more significant. For a companion with 
eh = 0.3, it reaches ~ 40% in the 20 - 30/im domain and, more 
importantly, it is enough to displace the total SED peak emis- 
sion from 15 to ~ 22/im. This effect is even more drastic in the 
ei, = 0.75 run, where the contribution of unstable grains balances 
that of < a cr j, particles for 25/mi< A < 60^/m. The increased 
importance of > a cr n grains for a more massive star is due to 
the fact that these grains are larger, because stiow ~ 3jum for an 
A3V star, and thus better emitters in the mid-IR than for an F5 
star. Furthermore, the star being brighter, the disc is hotter and 
the SED is shifted towards shorter wavelengths, where unstable 
grains are good emitters. For both et = 0.3 and et = 0.75 cases, 
the total SED's shift towards longer A is enough for it to peak at a 
wavelength slightly exceeding that of a ring of pa rticles located 
at a cr j, . In this case, an analysis similar to that of Trilling et al. 
(2007) would indeed infer the presence of dust in dynamically 
forbidden regions Q. 

Of course, our analysis is simplified. It neglects for instance 
the fact that grains beyond a c „, might not have a temperature 
corresponding to the radial distance at which they are observed. 
These grains are indeed on highly eccentric orbits covering a 
wide radial range and might not have the time to reach thermal 
equilibrium at a given distance r. This issue exceeds the scope 
of the present study, but our analysis has to the very least shown 
that grains in the dynamically unstable regions can significantly 
affect the disc's total SED, possibly making it look as if coming 
regions beyond 

4.1.2. Resolved observations: HR4796 

Our results are also applicable to resolved observations in scat- 
tered light, as they provide a tool to link the observed spatial 
structures to the potential influence of a companion star. One di- 
rect application are cases epitomized by HR4796, where a com- 
panion star, whose orbital parameters are unconstrained, is de- 
tected at a large projected distance from a circumprimary disc 
displaying pronounced features. As an example, HR4796B (at 
a projected distance ~510AU) has been invoke d as a potential 
source for HR4796A's disc sharp outer edge. [Augereau et alj 
(1999) postulated that if HR4796B's periastron is located at 
~ 200 AU, corresponding to et, ~ 0.4, then it could efficiently 
truncate the circumprimary disc at 70 AU. Howev er, these esti- 
mates were based on simulations by lHalll d 19971) of the early 
truncation of protoplanetary discs that neglected the effect of ra- 
diation pressure. 

We reinvestigate this scenario more thoroughly to test 
whether we can get a reasonable fit to the observed scat- 
tered light profile of the HR4796a disc with our more com- 
plete numerical model. We follow the procedure described in 



1.00 



- / / 



0.10 



0.01 




. Observed profile 

e,= 0.45 a 1 =350AU 
. e„=0.55 a„=455AU 
. e b =0.7 a B +760AU 
. e,=0.45 T = 5x1Cr' 



50 60 70 80 90 100 110 120 
Radial distance (AU) 

Fig. 9. HR4796A case. Steady state radial luminosity profile ob- 
tained for 4 different orbits of the companion and the observa- 
tionally derived estimate for the disc's optical depth r = 7.5 x 
10~ 3 , as well as for one alternative t = 5 x 10~ 4 value. The solid 
line i s the scattered light profile observed by Sch neider et al.l 
d2009l 



Sec l2.ll and first run a series of parent body simulations. We 
find that the minimum eccentricity for a perturber passing at 
r > 510AU to be able to truncate the parent body disc at 

~ 70 AU is eh — 0.45. This value is slighly higher than that 

i 1 j "~n 

of Augereau et al. ( 1999) because we assume the real mass ra- 

tio fx = 0.2 dJavawardhanaetal.l ll998) for t he HR4796 AB bi- 
nary instead of the fiducial n = 0.5 taken bv lHalll (1 19971) . Fig|9] 
shows the synthetic radial luminosity profile, in scattered light, 
obtained for 3 different binary configurations, ranging from et = 
0.45 to et = 0.7, and an optical depth t = 7.5 x 1 0~ 3 , an interme- 
diate v alue between the t 10~ 2 deduced from Au gereau et al.l 
(119991) and the t = 5 x 10 3 given by IWvatt et alj (fl999). 
Scattered light luminosity is derived assuming grey scattering 
of the grains. 

As can be clearly seen, the obtained luminosities in the 
r > a cr j t region are always much higher than the observed 
one, even for the lowest possible eccentricity value et = 0.45. 
However, a relatively satisfying fit can be obtained when lower- 
ing the disc's optical depth to t — 5 x 10~ 4 . This case might ap- 
pear to lack physical relevance as it requires an optical depth ap- 
proximately 15 times lower than the observed one. Nevetheless, 
this solution can not be directly ruled out because t is in our 
model mostly a way to parameterize the system's collisional ac- 
tivity. Even if there is usually a direct proportionality relation be- 
tween r and t co u, there exists some cases where this is not true for 
the relation between t and the collisional destruction timescale 
tDcoii- This might for instance be the case for dynamically cold 
discs, where impact velocities are too low to lead to shattering 
of the impactors, so that to C oll would be much higher than t co u 
deduced from t. Could HR4796 correspond to such a dynami- 
cally cold case? The answer to this question depends on the value 
of the proper eccentricity e p of particles within the ring, since 



7 Although iTrilling et al] i2007h used a simplified model, deriving 
grain characteristic temperatures under the assumption that they behave 
like black bodies. 



8 We leave the profiles in the r < «„,, region, especially inside the in- 
ner edge of the parent body ring at ~ 60 AU, out of our discussion. The 
sharp inner edge of the HR4796A disc must indeed have a specific ori- 
gin, unrelated to the companion, which is not the purpose of the present 
study. 
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dVcoii ~ SpVicep- Assuming a conservative threshold velocity for 
grain destruction of 50m.s _1 , this means that, at 70 AU from 
the primary, we need e p < 0.01. The only observational con- 
straint on ep is that the ring's width, 14 AU, gives an upper limit 
for e p equal to ~ 0.1. A e p < 0.01 solution is thus in principle 
possible, but it might not hold from a theoretical point of view. 
Debris discs are indeed systems in which the bulk of planetesi- 
mal accretion process should already be over and planets or plan- 
etary embryos are present and dynamically excite all remaining 
particles to eccentricities in the 0. 05 - 0.1 range (Artymowicz 
ll997tlThebault & Augereau II2007I) . As an example iWvatt et alJ 
(11999|) assumed (e p ) = 0.13 for their model of the HR4796A 
disc. As a consequence, the low collisional activity solution, 
even it cannot be fully dismissed, does not appear to be the most 
generic one. Furthermore, even if the disc is dynamically cold, 
then a sharp outer edge might be maintained without the need for 
an external perturber, as low e p val ues might naturally cre ate a 
depletion of small high-/? grains (see lThebault & Wul 2008). The 
companion is thus probably not a satisfying explanation for the 
HR4796A disc's radial profile. 

Our conclusion differs from that of Augereau e t~aT1 d 1999b 
because Hall (1997) did not consider the effect of radiation pres- 
sure, implicitly considering particles large enough to be only af- 
fected by gravity. This would correspond to the "parent body" 
disc of our study, for which we indeed find that it is truncated b y 
companion stars having similar characteristics as in 
But the present numerical exploration has shown that, even if 
HR4796B is eccentric enough to cut-off the circumprimary disc 
of parent bodies, it cannot do the same job for the small debris 
steadily produced by collisions amongst these parent bodies. 

More generally, these results might be useful for all systems 
where a bound, or possibly bound companion is detected far 
away from a circumprimary debris disc and is considered as a 
potential source for pronounced features in the disc, especially 
sharp edges, under the assumption that it is on a highly eccen- 
tric orbit allowing it to pass in the disc's vicinity. The present 
numerical exploration shows that the effect of such a compan- 
ion should in fact be limited, because companions on eccentric 
orbits are less efficient in preventing matter from populating dy- 
namically unstable region. 

4.2. Small grain depletion in dynamically stable regions 

In apparent contrast to the previous results showing the reduced 
effect of companion stars in the dynamically unstable regions, 
we find that binarity can also significantly affect a debris disc in 
the dynamically stable region, depleting it from a large fraction 
of its small grain population. This is once again due to radia- 
tion pressure: small particles produced at a < a cr u spend a large 
fraction of their eccentric orbits beyond a crit , where they can be 
removed by companion star perturbations. As a net result, small 
grains with sizes smaller than 3sthw (radiation pressure blow 
out size), will be underabundant in the collisionally active par- 
ent body region. This result is significant even if the parent body 
region lies well inside a„„, as shown in our inner ring case dis- 
played in Fig|7] 

Can this lack of very small grains have an observable coun- 
terpart? A priori, this should be the case for wavelengths com- 
parable to a given system's Sbiow Indeed, in the normal unper- 
turbed case, these discs' total cross section should be domi- 
nated by t he smallest bound grai ns, close to Sbiow (see for exam- 
ple [Thebault & Au gereaull2007l) . so that the disc should appear 
red for multi- wavelengths observations around A =* Sbiow This 



seem s for instance to be the case for HR4796A (ISchneider et al.l 
120091) . and might thus be an additional argument to rule out 
the binary sculpting scenario. Conversely, in the absence of 
Sbiow < 5 < 3 Sbiow grains, the disc should appear "grey" with 
no clear colour trend. 

Nevertheless, even if this effect is in principle qualitatively 
observable, these results should be taken with caution. Firstly 
because the size range for grain depletion is relatively limited 
(a factor 3), so that precise and well sampled (in wavelengths) 
observations should be necessary. In addition, different grain 
composition and porosity might greatly complicate this picture, 
changing the dynamical behaviour of a particle of a given size. 
Finally, there might be other mechanisms that can lead to a de- 
pl etion of small grains (a s for instance the dynamically cold case 
in lThebault &"^l2008h . 

5. Summary and Conclusions 

We have numerically explored the response of a circumprimary 
debris disc to the presence of a stellar companion. Our main re- 
sults can be summarized as follows: 

- A companion star is never able to fully truncate a a colli- 
sionally active circumprimary disc. Even if its perturbations 
place a critical distance a cr u to the primary beyond which or- 
bits are unstable, the r > a crlt regions are not empty. They are 
populated by small grains, steadily produced by collisions 
from parent bodies in the stable regions, which are placed on 
eccentric orbits by radiation pressure 

- The amount of material in the unstable regions depends on 
the balance between the rate at which small grains are pro- 
duced by collisions and the rate at which they are removed 
by the companion's perturbations. Massive discs with short 
collisional timescales are thus the most favourable cases, for 
which the system's total cross section can become dominated 
by the unstable grains. Conversely, for a fixed a crn , binaries 
with circular orbits are more efficient in removing r > a cr j, 
grains than binaries on highly eccentric orbits. 

- For massive, A-type primary stars, the contribution of the 
unstable grains can significantly affect the disc's SED. For 
dense discs and/or eccentric binaries, the global SED can 
peak at a wavelength longer than that corresponsping to 
grains located at a crU . 

- For discs that are resolved in scattered light and display sharp 
outer edges, our study shows that it is difficult to explain 
these outer truncations by a distant companion star on an ec- 
centric orbit. For the specific case of HR4796A, we show 
that the disc is probably too massive (i.e., collisionally ac- 
tive) and the detected companion on a too-eccentric orbit to 
explain the observed structure of the circumprimary disc. 

- Further putting into perspective the role of a e ,„ as an absolute 
frontier is the fact that the companion star is able to affect the 
grain population within the stable region. The r < a cra region 
is depleted from grains within a factor 2-3 of the blow out 
size by radiation pressure. However, this absence of small 
particles might not be easily observable. 

Acknowledgements. PT wishes to thank Mark Wyatt for enlightening discus- 
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